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SUMMARY 


A numerical procedure for the relaxation solution of the full steady Euler 
equations is described. By embedding the Euler system in a second-order surrogate 
system, central differencing may be used in subsonic regions while retaining matrix 
forms well suit'.'d to iterative solution procedures and convergence acceleration 
techniques. Hence, this method allows the development cf stable, fully-conservative 
differencing schemes for the solution of quite general inviscid flow problems. 
Results are presented for both subcritical and shocked, supercritical internal 
flows. Comparisons are made with a standard time-dependent solution algorithm. 


INTRODUCTION 


Techniques for the time-accurate solution of the unsteady Euler equations are 
well known and have relatively firm theoretical basis. The methodology for 
obtaining steady solutions to the Euler equations is in a more formative state and 
has undergone considerable evolution in the recent past. 

Given steady boundary conditions and assianing that a unique steady solution 
exists, one may solve the unsteady Euler equations in a non-tlme-accurate fashion by 
means of an algorithm with a stability bound in excess of the CFL condition or by 
using local time stepping. These constitute the simplest form of pseudo-unsteady 
solution procedure. Other oethods fix the total enthalpy at its steady-state value 
or use enthalpy damping, as proposed by Jameson, Schmidt and Turkel (1981), to 
accelerate convergence. Ni (1981) and Steger (1981) have developed multiple grid 
schemes for pseudo-unsteady solution of the time-dependent Euler equations. 
Approaches which modify the unsteady Euler equations by adding either time-dependent 
terms or time-dependent equations which enhance convergence have been developed by 
Essers (1980) and Viviand (1981). 

Attempts to deal directly with the steady first-order Euler system meet with 
immediate difficulties. Only centered difference operators will lead to numerical 
schoies which are simultaneously stable for the entire 8yst» in subsonic flow. 
However, such schemes lead to ill-conditioned matrices which defeat most Iterative 



solution procedures. 


Semi'slirect methods using Newton iteration to avoid this numerical difficulty 
have been developed by Rtzzi (1979) for the homoenthalpic Euler equations, by Rizzi 
and Skbllermo (1981) for the full potential equation written as a first-order 
system, and by Shubin, Stephens and Glaz (1981) for the quasi -one-dimensional Euler 
equations. The large matrices used by such methods may be an impediment when they 
are generalized to the full Euler system, or to higher dimensions. 

An alternative approach to solving the steady Euler equations involves either 
embedding them in a higher-orde’’ steady system or preconditioning their finite 
difference representations to enhance their compatibility with iterative procedures. 
Chattot, Guiu-Roux and Laminie (1981) used a variational approach to transform the 
first-order system representing the potential equation into an equivalent second- 
order system which was solved by a conjugate gradient algorithm. Lomax (1981), 
Desideri and Lomax (1981) and Lomax, Pulliam and Jespersen (1981) described a strat- 
egy for solving the Euler equations by preconditioning the finite difference equa- 
tions, choosing a stable relaxation procedure, and accelerating its convergence by a 
multiple grid technique. Preconditioned relaxation solutions were obtained for two- 
dimensional subsonic flow by Desideri and Lomax, while Lomax, Pulliam and Jespersen 
used a multiple grid procedure to solve a quasi-one-dimensional supersonic flow. 

In Johnson (1981) wo developed an approach whereby a first-order partial 
differential system Is embedded In a second-order surrogate system which may then be 
solved by means of the same sort of numerical techniques routinely used on the 
potential equation. We obtained results with the full Euler system for both 
supersonic and subcrltical, subsonic flow and with the transonic small disturbance 
equations for both subcrltical and supercritical flow. The embedding used with the 
Euler equations assumed the invertability of flux-vector Jacobian matrices and was 
thus not suitable for use in transonic flow computations. The present paper presents 
the details of an Improved surrogate equation technique which is capable of treating 
the full steady transonic Euler equations. 

SURROGATE EQUATION TECHNIQUE 

Given a first-order partial differential system, we embed this system in a 
second-order surrogate system, apply additional constraints to restrict the solution 
set of the surrogate, and solve the resulting paitial differential problem by means 
of a conventional iterative procedure. This method maintains the generality of the 
Euler equations, while allowing the use of the same sort of relaxation procedures 
developed for the efficient solution of second-order equations. 
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Consider a f irst-v)rdt>r system written in conservation law form, such as 

i (*)+^(>)] ,-0 

where q is an n-compunent vector and A and B are n x n matrices. We embed this 
system in a second-order surrogate of the form 

This system preserves the conservation law form of the o-iginal first-order system. 

The nature of the second-order partial differential operat->r is controlled by the 

choice of the matrices M and M. For example, the choice M <■ A^. N ■ B^ symmetrizes 

the coefficients of the terms of highest order and causes the sur&'ogate system to be 

T T 

elliptic, while the choice M ■ A , N ■ -B results In a non-elliptlc system which 
may be upwind differenced. Alternative choices for H and N are possible but will 
not be discussed here. The problem specification Is completed by requiring that, in 
addition to satisfying the original boundary conditions of the underlying 
first-order system, the solution to the surrogate system must also satisfy the 
first-order equations themselves at the boundaries. This is done to Insure 
uniqueness. Additionally, In the case where we employ an elliptic surrogate system 
to compute a supercritical flow, this boundary treatment allows the Introduction of 
dissipative terms for proper shock capture. 

Observe that, by switching the second-order operator from elliptic to 
hyperbolic type when the flow changes from subsonic to supersonic, it is possible to 
create a type-dependent differencing scheme for the surrogate system. Such a scheme 
could provide an alternative means for Insuring the correct shock capture and thus 
relax some of the constraints on the boundary treatment. While Initial computations 
indicate that this may be a viable approach, the results to be presented 
subsequent! V were obtained using the choice M • A , N ■ B everywhere in the domain. 
Notice that, in this case, the embedding operator is a close relative of the formal 
transpose of the Euler operator: 

A discrete representation of this operator has been independently proposed as a 
preconditioning operator by Desideri and Lomax. Because they operate on the finite 
difference equations, their approach is in several additional respects distinct from 
the one discussed here. For example, preconditioning the finite difference 
equations results in an effective non-compact differencing of the second-order 
system. Permutation of the resulting matrices is required to restore compact 
structure and reduce bandwidth. Furthermore, the surrogate equation technique 
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appears Co offer more flexibility in the treatment of boundary conditions than la 
available with the preconditioned finite difference equation approach. 

The second-order partial differential problem, being compatible with iterative 
techniques, may be solved by a variety of methods. For demonstration purposes, ve 
use f ully-conservative differencing together with the well-established successive 
line relaxation method. 


results 


We compute sober ideal and shocked, supercritical flows in a straight channel 
with a lOZ half-thick circular arc airfoil mounted on its lower wall. The 
second-order partial differencial problem is illustrated in Fig.l. As physical 
boundary conditions, we require that flow tangency be satisfied at solid walls, we 
specify total pressure, total temperature and flow angle at the inlet, and we 
specify the exit static pressure. We require that the first-order Euler equations 
be satisfied at each boundary to provide the additional boundary conditions 
necessary to completely pose the problem and to Insure the correct shock capture. 

As a standard of comparison for the accuracy of Che results presented here, we have 
recomputed all cases using the explicit MacCormack (1969) algorithm. 

The subcriclcal test case had an isentropic inlet Mach number of 0.5. Fig. 2 
shows the comparative upper and lower surface Mach number distributions. Isomach 
contours are plotted in Fig. 3. The supercritical test case was run at an inlet Mach 
number of 0.675, producing a shocked but unchoked flow. The comparative surface 
Mach number distributions and isomachs are shown in Figs. 4 and 5, respectively. The 
sonic line is dashed in the isomach plots. Comparison of the results of the 
surrogate equation algorithm with those of the MacCormack algorithm is encouraging. 
Minor discrepancies may be attributed, in part, to differences in the truncation 
error of the two algorithms or to the lack of complete annihilation of low frequency 
error modes. 

Representative convergence histories are shown in Fig. 6. The surrogate 
equation algorithm, using successive line overrelaxation on the second-order system, 
converges more rapidly than the MacCormack algorithm, using local time stepping at 
0.9 of the local CFL limit. Since the residuals are defined quite differently for 
the tvo methods, the error measure used for this comparison is the correction to the 
vector of conservation variables. Consequently, the relative convergence rates are 
more significant than the indicated levels of error. Having demonstrated the 
capability of the surrogate equation technique to solve the full Euler equations by 
relaxation, it should prove relatively straightforward to further accelerate 
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coiwiTgi-nco. 


W)iiU' all reBultd prcBimtrd here are two-dimensional, the extension of the 
aurroKate equation technique to three dimenslona presents no essential difficulties. 
Furthermore, the use of this technique as an Invlscld component of a aonal procedure 
for the iterative solution c( the steady Navler-Stokes equations appears feasible. 


C ONCLUSIONS 

We nmy obtain a solution to the full steady transonic Euler equations by using 
a surrogate second-order system together with the original Euler physical boundary 
conditions and additional constraints obtained from the first-order Euler system. 

Tills surrogate equation technique provides a means for formulating problems 
involving the full steady Euler equations in such a way as to allow the use of 
stable, fully-conservatlve differencing and relaxation solution procedures. Hence, 
we may solve cither Irrotatlonal or rotational flow problems across the entire 
spectrum of subsonic, transonic and supersonic conditions without resort either to 
derived dependent variables, semi-direct methods, or to an unsteady formulation. 

^bedding the Euler equations in a second-order system allows the application 
of the many convergence acceleration techniques which have been developed for other 
second-order systems. Thus, the surrogate equation technique provides an 
opportunity for the construction of fast and efficient numerical procedures for the 
solution of the full steady Euler equations. 
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Fig. 1 Partial Differential Problem 



Axial Coordinate Axial Coordinate 


Surrogate Equation Algorithm MacCormack Algorithm 

Fig. 2 Subcritical Mach Number Distributions 
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Fig. 3 Subcritical laomach Contours 
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Fig. 5 Supercritical laomach Contoura 
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Fig. 6 Convergence Histories 





